Onset temperature of Bose-Einstein condensation in incommensurate solid 4 He 
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Abstract 

The temperature dependence of the one-body density matrix in 4 He crystals presenting vacancies is com- 
puted with Path Integral Monte Carlo. The main purpose of this study is to estimate the onset temperature 
To of Bose-Einstein condensation in these systems. We see that To depends on the vacancy concentration X v 

2/3 

of the simulated system, but not following the law To ~ X v obtained assuming non-interacting vacancies. 
For the lowest X v we have studied, that is X v — 1/256, we get To = 0.15 ±0.05 K, close to the temperatures 
at which a finite fraction of non-classical rotational inertia is experimentally observed. Below To, vacancies 
do not act as classical point defects becoming completely delocalized entities. 

PACS numbers: 67.80.-s, 61.72.Bb 
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The debate about the supersolid state of matter, i.e., a phase where crystalline order coex- 
ists with superfluidity, has gained a great interest among the scientific community after the first 
observation of non classical rotational inertia (NCRI) in torsional oscillators containing solid he- 
lium [H, Q|. Although several experiments have confirmed the appearance of a phase transition in 
solid 4 He at temperatures T c ~ 60 — 100 mK 0,0], we are still far from a complete description 
of this phenomenon because of controversial experimental results. For instance, the values of the 
superfluid density p s /p reported so far can vary more than one order of magnitude according to 
experimental conditions such as the way in which the crystal is prepared, its subsequent annealing 
or the 3 He concentration 0-0]- These discrepancies suggest that the quality of the solid sample 
plays a very important role in these experiments and make fundamental a study of crystalline 
defects in quantum crystals. 

First theoretical studies, indeed, suggested that a possible supersolid behavior can be explained 
assuming the presence, in the ground state of quantum crystals, of delocalized vacancies which may 
undergo Bose-Einstein condensation (BEC) at low temperature 0]. Nevertheless, these early works 
were based on simplified models, so that it was not possible to draw specific predictions for solid 4 He. 
More recently, microscopic methods have been extensively used to provide a reliable description 
of the supersolid state but, so far, they have not been able to reproduce all the experimental 
findings. Path integral Monte Carlo (PIMC) simulations have shown that a commensurate perfect 
crystal does not exhibit superfluidity 0-11], but a non-zero condensate fraction has been observed 
in crystals with a finite vacancy concentration at zero temperature 121 ]. 

The possibility for solid 4 He to present vacancies in its ground state seems to be hindered by the 
energetic cost of these defects. According to several Quantum Monte Carlo results, the vacancy 



formation ener gy is estimated to be of the order of 10 K 13l4l71]. in agreement with experimental 
measurements 18y. Nevertheless, the high derealization of the vacancies in solid 4 He at tempera- 



tures close to zero prevents an interpretation of these defects in terms of a classical theory involving 
an activation energy and a configurational entropy for their creation Hi- 21]. Furthermore, experi- 
mental data cannot rule out the possibility of a zero-point vacancy concentration below 0.4% [22]. 
It has also to be noticed that formation energy considerations do not exclude the possibility of 
vacancies introduced through the experimental conditions, for example during the crystal growth. 
The spatial correlation between vacancies has been calculated in order to understand if a gas of 
defects can be metastable in solid 4 He. The results show an attractive correlation between vacan- 
cies at short distance, but they cannot conclude if they form bound states and aggregate in large 
clusters which eventually would phase separate 

ass 

In this work, we calculate by means of the Path Integral Monte Carlo (PIMC) method, the 
one-body density matrix pi(v, r') in solid 4 He samples presenting a finite vacancy concentration, 
focusing especially on its temperature dependence. In the study of the BEC properties of quantum 
systems, pi(r, r') is a fundamental quantity, the condensate fraction no being its asymptotic limit 
for large [r — r'| values. Our main purpose is to estimate the onset temperature of BEC To and 
to compare it with the experimental measurements. We start simulating an hep crystal with a 
vacancy concentration X v = 1/180, trying also to give a qualitative picture of the derealization 
of the vacancies and the appearance of BEC. Finally, we study the dependence of Tq on X v to 
guess which would be the vacancy concentration needed to have BEC appearing in the range of 
temperatures T c ~ 60 — 100 mK at which NCRI is experimentally observed. 

PIMC provides a fundamental approach in the study of the thermodynamic properties of 
strongly interacting quantum systems at finite temperature [26[. In this method, the partition 
function Z is rewritten making use of the convolution property of the thermal density matrix 
G(R' , R; P) = (R'\e~P H \R) (where j3 = {ksT)~ l is the inverse of the temperature and H is the 
Hamiltonian of the system), which is known only for small ft. This procedure is equivalent to map- 
ping the quantum many-body system at finite temperature onto a classical system of closed ring 



polymers. Increasing the number of convolution terms used to rewrite Z, which corresponds to the 
number of beads composing each classic polymer, one is able to reduce the systematic error due 
to the approximation for G and therefore to recover "exactly" the thermal equilibrium properties 
of the system. A good approximation for the propagator G is fundamental in order to reduce the 
complexity of the calculation and ergodicity issues. Using the Chin approximation 271. l28l| . we are 
able to obtain an accurate estimation of the relevant physical quantities with reasonable numeric 
effort even in the low temperature regime, where the simulation becomes harder due to the large 
zero-point motion of particles. Chin approximation for the action is accurate to fourth order in 
the imaginary-time step but a real sixth-order behavior can be achieved by adjusting properly the 
two parameters entering in it. Similar accuracies can be achieved using other high-order proposals 
for the action 29T-31]. 



An additional problem we have to deal with when approaching the low temperature limit with 
PIMC simulations arises from the indistinguishable nature of 4 He atoms. Since we study a bosonic 
system, the symmetry of Z can be recovered via the direct sampling of permutations between the 
ring polymers. To this purpose, we have used the Worm Algorithm [32]. This algorithm allows 
PIMC simulations with a very efficient sampling of the exchanges between bosons. Furthermore, it 
is able to give an estimation of the normalization factor of pi(ri, avoiding thus the systematical 
uncertainties which can be introduced by a posteriori normalization factor. 

In order to calculate pi(r) = Pi(|r — r'|) in a crystal with vacancy concentration X v = 1/180, 
we have carried out simulations of N = 179 4 He atoms, interacting through an accurate Aziz 
pair potential [jjij, in an almost cubic simulation box matching the periodicity of an hep lattice 
made up of N s = 180 sites at a density p = 0.0294 A -3 . We apply periodic boundary conditions 
to the simulation box in order to simulate the infinite dimensions of the bulk system. In Fig. 
[H we show the results for pi(r) at different temperatures and we compare them with the zero 
temperature estimations of p\ for the same system and for a perfect hep crystal, obtained with 
the Path Integral Ground State method in Ref. [3J]. We notice that, at temperatures T > 0.75 
K, p\{r) computed in an incommensurate crystal, even though is not compatible with p\ for the 
perfect crystal, presents a similar exponential decay at large r. At lower temperatures, the decay 
of pi(r) is smoother and, for temperatures T < Tq = 0.2 K, p\ presents a non-zero asymptote 
at large r, which indicates the presence of BEC inside the system. This To can be considered a 
first estimate of the onset temperature of supersolidity in the simulated system. An analysis of 
the finite size effects would be needed to get a more precise estimation of the critical temperature 
of the supersolidity transition. Nevertheless, the simulation of bigger systems with exactly the 
same vacancy concentration requires a huge computational effort that would make the calculations 
impracticable. 

In order to give a more qualitative description of the appearance of BEC in incommensurate 4 He 
solids, we visualize typical configurations of the system during the simulation. In Fig. [21 we plot 
two-dimensional projections of the positions of the quantum particles (represented by polymers 
in PIMC) lying in a basal plane of the incommensurate hep crystal at different temperature. At 
T = 1 K, 4 He atoms tend to stay localized around their equilibrium positions. Also the vacancies 
are localized and can be easily detected inside the lattice. This explains the fact that, at that 
temperature, the presence of vacancies does not affect noticeably the overall behavior of p\ which, 
for the incommensurate crystal, is similar to the one of the perfect crystal. At T = 0.5 K, the 
effects of the derealization of the 4 He atoms can be seen with the appearance of some polymers 
which are spread on two different lattice points. In the space configurations at this temperature, 
the acceptance rate of the exchange between the polymers is higher than in the configurations 
at larger temperature, but it is still too low to allow the appearance of long permutation cycles, 
which are necessary to see BEC. At T = 0.2 K, the large zero-point motion of the He atoms 
makes the vacancy delocalized and undetectable inside the crystal, which looks like a commensurate 
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Figure 1: (Color online) The one-body density matrix pi(r) for an hep crystal with vacancy concentration 
X v = 1/180 at density p = 0.0294 A" 3 at different temperatures: T = 1 K (red triangles down), T = 0.75 K 
(blue triangles up), T = 0.5 K (green diamonds), T = 0.25 K (yellow squares) and T — 0.2 K (purple circles). 
The dotted and dashed lines represent p\(r) at zero temperature respectively for the commensurate (X v = 0) 
and incommensurate crystal (X v — 1/180) at the same density, taken from Ref. 34 1. 



system. Since the number of lattice sites is different from the number of particles, this means that 
different polymers may superpose over the same lattice site: this occurrence strongly enhances the 
possibility for the atoms to permutate and allows the creation of long permutation cycles which 
close on periodic boundary conditions. The appearance of configurations presenting a non zero 
winding number, as the one shown in Fig. [3j indicates that the simulated crystals below To = 0.2 
K support superfluidity. However, it is not possible to give a reliable estimation for the superfluid 
density p s /p in these systems, since the smallest value for p s /p computable with the winding 
number estimator is of the order of 1%, that is of the same order of the value expected from the 
experimental measurements. 

In order to study how the vacancy concentration in quantum solids affect the onset temperature 
of BEC, we have computed the one-body density matrix also for fee 4 He crystals with X v = 1/108, 
X v = 1/128, and X v = 1/256. It is worth noticing that the result for X v = 1/128 has been 
obtained in a simulation with two vacancies in a lattice of 256 site points and the results for both 
the condensate fraction and onset temperature follow the same X v dependence that the single 
vacancy cases. In Table [U we show the onset temperature of BEC To and the condensate fraction 
no at low temperature obtained with PIMC in the three crystals we have studied. We notice that, 
for the lowest X v , we get To = 0.15±0.05, which is close to the temperatures at which supersolidity 
has been experimentally observed. 

In Fig. [H we plot our results for To as a function of X v . Our results for To do not follow the 
law To ~ X^ 3 , obtained from a description of solid 4 He in terms of a rarefied Gross-Pitaevskii 
superfluid gas of vacancies, as proposed by Anderson in Ref. [2lJ. This seems to suggest that, 



at least in the range of X v we have been able to study, the correlations between vacancies have 
an important effect on Tq and the system cannot be described within a mean-field approach. 
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Figure 2: (Color online) Two-dimensional projection of basal planes of the incommensurate hep crystal at 
different temperatures, represented according to the PIMC isomorphism of the classical polymers. At T = 1 
K (higher panel) the vacancy is localized and indicated by the red circle. At T — 0.5 K (middle panel), the 
vacancy begins to delocalize: the red ellipse indicate a quantum particle delocalized over two different lattice 
sites. At T = 0.2 K (lower panel), the vacancy is completely delocalized and cannot be easily detected. 
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Figure 3: (Color online) Two-dimensional projection of two consecutive basal planes of the incommensurate 
hep crystal at T = 0.15 K. The different colors distinguish the two different planes. The thick solid line 
represent a long permutation cycle between the 4 He atoms presenting a non-zero winding number. 



N v 


N s 


X v 


To (K) 


n Q 


1 


108 


1/108 


0.50 ±0.10 


(1.81 ±0.14)10~ 3 


2 


256 


1/128 


0.40 ± 0.075 


(1.09±0.13)10- 3 


1 


180 


1/180 


0.20 ±0.05 


(9.0±0.8)10- 4 


f 


256 


1/256 


0.15 ±0.05 


(7.2±0.8)10- 4 



Table I: The onset temperature of BEC To and the condensate fraction no at low temperature as a function 
of the vacancy concentration X v — N v /N s , N v and N s being the number of vacancies and of lattice sites, 
respectively. 



Nonetheless, our qualitative description of 4 He crystals supports the hypothesis [2l| according to 
which it is not reasonable to regard vacancies in quantum solids as strictly local entities. 

In an attempt to estimate which should be the vacancy concentration in 4 He crystals needed to 
have BEC appearing at the temperature T c measured experimentally for the supersolid transition, 
we have plotted in Fig. U] a power function trying to fit the PIMC results. According to this 
empirical law, 4 He crystals with a vacancy concentration X v ~ 2 — 3 x 10~ 3 would have an onset 
temperature To in agreement with the experimental values T c . This result for X v is in good 
agreement with the equilibrium vacancy concentration in solid 4 He at zero temperature obtained 
variationally with the shadow wave function [24f ] . 

In conclusion, we have shown that the onset temperature To of BEC in 4 He crystals presenting 
vacancies, calculated using the PIMC method, is comparable with the experimental measurements 
of the supersolid transition temperature when the concentration of vacancies is small enough (X v ~ 
2 — 3 x 10 -3 ). PIMC simulations also show clearly that when this onset temperature is reached the 
vacancies become completely delocalized objects, as hypothesized in the past [19|, l2l| and never 
microscopically observed so far. 

This work was partially supported by DGI (Spain) under Grant No. FIS2008-04403 and Gen- 
eralitat de Catalunya under Grant No. 2009-SGR1003. 
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Figure 4: (Color online) The onset temperature To of BEC in incommensurate solid 4 He as a function of 
the vacancy concentration X v : the blue squares represent the results obtained with the PIMC method; the 
red line is a fit to the data with a power law To = AX® with optimal values A = 789 K and B = 1.57; the 
grey band indicates the temperature at which the NCRI appears experimentally. 
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